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Abstract. Declarative spatial reasoning denotes the ability to (declar- 
atively) specify and solve real-world problems related to geometric and 
qualitative spatial representation and reasoning within standard knowl¬ 
edge representation and reasoning (KR) based methods (e.g., logic pro¬ 
gramming and derivatives). One approach for encoding the semantics of 
spatial relations within a declarative programming framework is by sys¬ 
tems of polynomial constraints. However, solving such constraints is com¬ 
putationally intractable in general (i.e. the theory of real-closed fields). 

We present a new algorithm, implemented within the declarative spatial 
reasoning system CLP(QS), that drastically improves the performance 
of deciding the consistency of spatial constraint graphs over conventional 
polynomial encodings. We develop pruning strategies founded on spatial 
symmetries that form equivalence classes (based on affine transforma¬ 
tions) at the qualitative spatial level. Moreover, pruning strategies are 
themselves formalised as knowledge about the properties of space and 
spatial symmetries. We evaluate our algorithm using a range of bench¬ 
marks in the class of contact problems, and proofs in mereology and 
geometry. The empirical results show that GLP(QS) with knowledge- 
based spatial pruning outperforms conventional polynomial encodings 
by orders of magnitude, and can thus be applied to problems that are 
otherwise unsolvable in practice. 

Keywords: Declarative Spatial Reasoning, Geometric Reasoning, Logic 
Programming, Knowledge Representation and Reasoning 


1 Introduction 


Knowledge representation and reasoning (KR) about space may be formally in¬ 
terpreted within diverse frameworks such as: (a) analytically founded geometric 
reasoning & constructive (solid) geometry [21, 27, 29]; (b) relational algebraic se¬ 
mantics of ‘qualitative spatial calculi’ [24]; and (c) by axiomatically constructed 
formal systems of mereotopology and mereogeometry [1[. Independent of for¬ 
mal semantics, commonsense spatio-linguistic abstractions offer a human-centred 
and cognitively adequate mechanism for logic-based automated reasoning about 
spatio-temporal information [5]. 

[> Declarative Spatial Reasoning In the recent years, declarative spatial 
reasoning has been developed as a high-level commonsense spatial reasoning 


2 


Carl Schultz and Mehul Bhatt 


paradigm aimed at (declaratively) specifying and solving real-world problems 
related to geometric and qualitative spatial representation and reasoning [4]. A 
particular manifestation of this paradigm is the constraint logic programming 
based CLP(QS) spatial reasoning system [4, 33, 34] (Section 2). 

[> Relational Algebraic Qualitative Spatial Reasoning The state of the 
art in qualitative spatial reasoning using relational algebraic methods [24] has 
resulted in prototypical algorithms and black-box systems that do not integrate 
with KR languages, such as those dealing with semantics and conceptual knowl¬ 
edge necessary for handling background knowledge, action & change, relational 
learning, rule-based systems etc. Furthermore, relation algebraic qualitative spa¬ 
tial reasoning (e.g. LR ]25]), while efficient, is incomplete in general ]22, 23, 24]. 
Alternatively, constraint logic programming based systems such as CLP(QS) ]4] 
and others (see ]9, 10, 18, 20, 28, 29]) adopt an analytic geometry approach where 
spatial relations are encoded as systems of polynomial constraints;® while these 
methods are sound and complete (see Section 2.2), they have prohibitive com¬ 
putational complexity, 0(c)i^) in the number of polynomial variables n, meaning 
that even relatively simple problems are not solved within a practical amount 
of time via “naive” or direct encodings as polynomial constraints, i.e. encodings 
that lack common-sense knowledge about spatial objects and relations. On the 
other hand, highly efficient and specialised geometric theorem provers (e.g. ]12]) 
and geometric constraint solvers (e.g. ]17, 27]) exist. However, these provers ex¬ 
hibit highly specialised and restricted spatial languages® and lack (a) the direct 
integration with more general AI methods and (b) the capacity for incorpo¬ 
rating modular common-sense rules about space in an extensible domain- and 
context-specific manner. 

The aims and contributions of the research presented in this paper are two-fold: 


1. to further develop a KR-centered declarative spatial reasoning paradigm 
such that spatial reasoning capabilities are available and accessible within 
AI programs and applications areas, and may be seamlessly integrated with 
other AI methods dealing with representation, reasoning, and learning about 
non-spatial aspects 


^ Incompleteness refers to the inability of a spatial reasoning method to determine 
whether a given network of qualitative spatial constraints is consistent or inconsistent 
in general. Relation-algebraic spatial reasoning (i.e. using algebraic closure based 
on weak composition) has been shown to be incomplete for a number of spatial 
languages and cannot guarantee consistency in general, e.g. relative directions ]23] 
and containment relations between linearly ordered intervals [22], Theorem 5.9. 

® We emphasise that this analytic geometry approach that we also adopt is not qual¬ 
itative spatial reasoning in the relation algebraic sense; the foundations are similar 
(i.e. employing a finite language of spatial relations that are interpreted as infinite 
sets of configurations, determining consistency in the complete absence of numeric 
information, and so on) but the methods for determining consistency etc. come from 
different branches of spatial reasoning. 

® Standard geometric constraint languages of approaches including ]12, 17, 27] consist 
of points, lines, circles, ellipses, and coincidence, tangency, perpendicularity, paral¬ 
lelism, and numerical dimension constraints; note the absence of e.g. mereotopology 
and “common-sense” relative orientation relations ]35]. 
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2. to demonstrate that in spite of high computational complexity in a general 
and domain-independent case, the power of analytic geometric —in par¬ 
ticular polynomial systems for encoding the semantics of spatial relations- 
can be exploited by systematically utilising commonsense knowledge about 
spatial object and relationships at the qualitative level. 


We present a new algorithm that drastically improves analytic spatial reasoning 
performance within KR-based declarative spatial reasoning approaches by iden¬ 
tifying and pruning spatial symmetries that form equivalence classes (based on 
affine transformations) at the qualitative spatial level. By exploiting symmetries 
our approach utilises powerful underlying, but computationally expensive, poly¬ 
nomial solvers in a significantly more effective manner. Our algorithm is simple 
to implement, and enables spatial reasoners to solve problems that are otherwise 
unsolvable using analytic or relation algebraic methods. We emphasise that our 
approach is independent of any particular polynomial constraint solver; it can be 
similarly applied over a range of solvers such as CLP(R), SMTs, and specialised 
geometric constraint solvers that have been integrated into a KR framework. 

In addition to AI / commonsense reasoning applications areas such as design, 
GIS, vision, robotics [3, 5, 6, 7], we also address application into automating 
support for proving the validity of theorems in mereotopology, orientation, shape, 
etc. (e.g. [8, 36]). Building on such foundational capabilities, another outreach is 
in the area of computer-aided learning systems in mathematics (e.g. at a high- 
school level). For instance, consider Proposition 9, Book I of Euclid’s Elements, 
where the task is to bisect an angle using only an unmarked ruler and collapsable 
compass. Once a student has developed what they believe to be a constructive 
proof, they can employ declarative spatial reasoners to formally verify that their 
construction applies to all possible inputs (i.e. all possible angles) and manipulate 
an interactive sketch that maintains the specified spatial relations (i.e. dynamic 
geometry [17]). A further area of interest is verifying the entries of composition 
tables that are used in relation algebraic qualitative spatial reasoning ]30] : given 
spatial objects a,b,c G U, composition “look up” tables are indexed by pairs of 
(base) relations Ri^^, and return disjunctions of possible (base) relations 
For each entry, the task is to prove 3a,b,c G U[Ri^^ A i?2bc ^ -^3ac) 
only those base relations i ?3 in the entry’s disjunction. 


2 Declarative Spatial Reasoning with CLP(QS) 


Declarative spatial reasoning denotes the ability of declarative programming 
frameworks in AI to handle spatial objects and the spatial relationships amongst 
them as native entities, e.g., as is possible with concrete domains of Integers, 
Reals and Inequality relationships. The objective is to enable points, oriented 
points, directed line-segments, regions, and topological and orientation relation¬ 
ships amongst them as first-class entities within declarative frameworks in AI 

[4]. 
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2.1 Examples of Declarative Spatial Reasoning with CLP(QS) 


With a focus on spatial question-answering, the CLP(QS) spatial reasoning sys¬ 
tem [4, 33, 34] provides a practical manifestation of certain aspects of the declar¬ 
ative spatial reasoning paradigm in the context of constraint logic programming 
(CLP). ^ CLP(QS) utilises a high-level language of spatial relations and com- 
monsense knowledge about how various spatial domains behave. Such relations 
describe sets of object configurations, i.e. qualitative spatial relations such as 
coincident, left of, or partially overlapping. Through this deep integration of 
spatial reasoning with KR-based frameworks, the long-term research agenda 
is to seamlessly provide spatial reasoning in other AI tasks such as planning, 
non-monotonic reasoning, and ontological reasoning [4]. What follows is a brief 
illustration of the spatial Q/A capability supported by CLP(QS). 

EXAMPLE A. Massachusetts Comprehensive Assessment System (MCAS) 

Grade 3 Mathematics (2009), Question 12. Put a square and two right- 
angled triangles together to make a reetangle. (1) Put the shapes Ti,T 2 ,S illus¬ 
trated in Figure 1(d) together to make a rectangle. (2) Put the shapes Ti,T 2 ,S 
in Figure 1(d) together to make a quadrilateral that is not a rectangle. 

CLP(QS) represents right-angle triangles as illustrated in Figure 1(b). Figures 
1(a) and 1(c) present the CLP(QS) solutions. 

Grade 3 Mathematics (2013), Question 17. (1) Flow many copies of Ti 
illustrated in Figure 1(d) are needed to eompletely fill the region R illustrated in 
Figure 2(a) without any of them overlapping? 

As presented in Figure 2, CLP(QS) solves both the geometric definition and a 
variation where the dimensions of the rectangle and triangles are not given. 

EXAMPLE B. Qualitative Spatial Reasoning with Complete Unknowns 

In this example CLP(QS) reasons about spatial objects based solely on given 
qualitative spatial relations, i.e. without any geometric information. 

Define three cubes A, B, C. Put B inside A, and make B disconnected from C. 
What spatial relations can possibly hold between A and C ? 

CLP(QS) determines that A must be disconnected from C and provides the 
inferred corresponding geometric constraints, as illustrated in Figure 3. 


2.2 Analytical Geometry Foundations for Declarative Spatial 
Reasoning 

Analytic geometry methods parameterise classes of objects and encode spatial 
relations as systems of polynomial equations and inequalities [12]. For example, 
we can define a sphere as having a 3D centroid point {x,y,z) and a radius r, 
where x,y,z,r are reals. Two spheres si,S 2 externally connect or touch if 

(Xs, - Xs^f + (ysi - + {Zs^ - Zs^f = {vs, + (1) 

^ Spatial Reasoning (CLP(QS)). www.spatial-reasoning.com 
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?- T1 = right_triangle(_,_,W,W) , 

I T2 = right_triang‘le(_,_,W,W), 

I S = rectangle(_,W,W), 

I Solution = rectangle(_,_,_), 

I topology(rcc8(eq), Solution, tile_union(IT1,T2,S])), 

I ground_object(Solution). 

T1 = right_triangle(orientation(0) , point(0, 0), 5, 5), 

W = 5, 

T2 = right_triangle(orientation(180) , point(5, 5), 5, 5), 

S = rectangle(point(5, 0), 5, 5), 

Solution = rectangle(point(0, 0), IG, 5) . 

(a) CLP(QS) solution to arranging T\,T2, S to 
form rectangles. 

?- T1 = right_triangle(_,_,W,W) , 

I T2 = right_triangle(_,_,W,W) , 

I S = rectangle(_,W,W) , 

I topology ( rcc8 (eq), Solution, tile_union{ [T1,T2,S] )), 

I quadrilateral(Solution) , 

I not(Solution = rectangle(_,_,_) ), 

I ground_object(Solution). 

T1 = right_triangle{orientation(90), point(5, 0), 5, 5), 

W = 5, 

T2 = right_triangle(orientation(270) , point{10, 5), 5, 5), 

S = rectangle(point{5, 0), 5, 5), 

Solution = parallelogram(orientation{G), point(O,0), 10, 5); 

T1 = right_triangle{orientation(90) , point(5, 0), 5, 5), 

W = 5, 

T2 = right_triangle{orientation(0), point(10, 0), 5, 5), 

S = rectangle(point(5, G), 5, 5), 

Solution = trapezoid(orientation (0), point(0, G), 15, 5, 5); 
false. 

(c) CLP(QS) solution to arranging Ti, T2, S to 
form non-rectangular quadrilaterals. 


Fig. 1. Using CLP(QS) to solve MCAS 


Tw 


h.I 



orientation = 0° orientation = 90° 

(b) Representing right-angle triangle 
T in CLP(QS). 



shapes to arrange 



trapezoid 


(d) Shapes Ti,T2,S 
to be arranged and 
CLP(QS) solutions. 

3 Mathematics Test questions (2009). 


region R tiling solution 



?- between(l,8,TriangleCount), 

I right_tri3ngle_list(TriangleCount,[Tl|TRest]), 

I Tl = right_triangle(_,_,1,1), 

I size(equal, Tl, TRest), 

I topology(rcc8(eq), rectangle(point(G,0),3,1), 

I tile_union([Tl|TRest])). 


TriangleCount = 6 


?- between( 1,8, TriangleCount) , 

I right_triangle_list (TriangleCount , Triangles) , 

I topology( rcc8(eq) , rectangle(_,_,_) , 

I tile_union (Triangles )). 

TriangleCount = 2;... 

TriangleCount = 4;... 

TriangleCount = 6;... 

TriangleCount = 8; 
false. 


(a) CLP(QS) solution to tiling a rectangular 
region with triangle Ti. 


(b) When no geometric information is given, 
CLP(QS) determines that the number of right- 
angle triangles must be even. 


Fig. 2. Using CLP(QS) to solve MCAS Grade 3 Mathematics Test questions (2013). 
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?- A=cube(_, _) , B=cube{_, _), C=cube(_, _), 

I topology(rcc8(ntpp), A, B), 

I topology(rcc8(dc), B, C), 

I topology(Relation, A,C), 

I constraints([(A,’A'),(B,'B'),{C,'C’)],Cons), 

I pretty_print_constraints(Cons). 


CLP(QS) Constraints: 


(1) 

B.length+B.X'C.x< -Q.O 

(2) 

A.z-B.z>G.O 

(3) 

A. length- B .length+A. z-B. z<0 .G 

(4) 

A.y-B.y>G.O 

(5) 

A. length-B.length+A.y-B.y<0.G 

(6) 

A.x-B.x>G.G 

(7) 

A. length•B.length+A.x-B.x<0.G 

(8) 

A.length>0.G 

(9) 

C.length>0.0 


Relation = rcc8(dc), 





(1) B^ + B,< Cx 

(2) > B^ 

(3) + A, < B^ + Bi 

{4) A^ > By 

(5) + Ai < By + B; 

(6) A:, > 

{7) A^+Ai < B.^ + B/ 

(8) A, > 0 

(9) Q > 0 


B is left of C 

the front face of A is behind 
the front face of B 
the back face of A is in front 
of the back face of B 
the base of A is above the 
base of B 

the top of A is below the 
top of B 

the left side of A is right of 
the left side of B 
the right side of A is left of 
the right side of B 
A side length is positive 
C side length is positive 


(a) CLP(QS) solution and its corresponding (b) Geometric inequalities provided by 
inferred CLP(QS) that correspond to the qualitative 

geometric constraints. constraints. 


Fig. 3. Spatial reasoning about cubes A,B,C with complete geometric unknowns. 

If the system of polynomial constraints is satisfiable then the spatial constraint 
graph is consistent. Specifically, the system of polynomial (in) equalities over vari¬ 
ables X is satisfiable if there exists a real number assignment for each x G X such 
that the (in)equalities are true. Partial geometric information (i.e. a combina¬ 
tion of numerical and qualitative spatial information) is utilised by assigning the 
given real numerical values to the corresponding object parameters. Thus, we can 
integrate spatial reasoning and logic programming using Constraint Logic Pro¬ 
gramming (CLP) [19]; this system is called CLP over qualitative spatial domains. 
CLP(Q5), provides a suitable framework for expressing and proving first-order 
spatial theorems. 

Cylindrical Algebraic Decomposition (CAD) [13] is a prominent sound and com¬ 
plete algorithm for deciding satisfiability of a general system of polynomial con¬ 
straints over reals and has time complexity 0(cf ) in the number of free vari¬ 
ables [2[. Thus, a key focus within analytic spatial reasoning has been methods 
for managing this inherent intractibility.® More efficient refinements of the orig¬ 
inal CAD algorithm include partial CAD [14]. Symbolic methods for solving 
systems of multivariate equations include the Grobner basis method [11] and 
Wu’s characteristic set method [40] . In the QUAD-CLP(R) system, the authors 
improve solving performance by using linear approximations of quadratic con¬ 
straints and by identifying geometric equivalence classes [28] . Ratschan employs 
pruning methods, also at the polynomial level, in the rsolve system [31, 32]. 


® Important factors in determining the applicability of various analytic approaches are 
the degree of the polynomials (particularly the distinction between linear and non¬ 
linear) and whether both equality and inequalities are permitted in the constraints. 
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Constructive and iterative (i.e. Newton and Quasi-Newton iteration) methods 
solve spatial reasoning problems by “building” a solution, i.e. by finding a con¬ 
figuration that satisfies the given constraints [27]. If a solution is found, then the 
solution itself is the proof that the system is consistent - but what if a solution 
is not found within a given time frame? In general these methods are incomplete 
for spatial reasoning problems encoded as nonlinear equations and inequalities 
of arbitrary degree.® 

3 Spatial Symmetries 

Information about objects and their spatial relations is formally expressed as 
a constraint graph G = {N,E), where the nodes N of the graph are spatial 
objects and the edges between nodes specify the relations between the objects. 
Objects belong to a domain, e.g. points, lines, squares, and circles in 2D Eu¬ 
clidean space, and cuboids, vertically-extruded polygons, spheres, and cylinders 
in 3D Euclidean space. We denote the object domain of node i as Ui (spatial 
domains are typically infinite). A node may refer to a partially ground, or com¬ 
pletely geometrically ground object, such that Ui can be a proper subset of the 
full domain of that object type. Each element i' € Ui is called an instance of 
that object domain. A configuration of objects is a set of instances {i'l,... ,i'^} 
of nodes ii,... ,in respectively. 

A binary relation Rij between nodes i,j distinguishes a set of relative configura¬ 
tions of i,j\ relation R is said to hold for those configurations, Rij Q UiX Uj. In 
general, an n-ary relation for n > 1 distinguishes a set of configurations between 
n objects: U Ui^ x • • • x Ui^. 

An edge between nodes i,j is assigned a logical formula over relation symbols 
Ri ,..., Rjn and logical operators V, A, Given an interpretation i',j', the for¬ 
mula for edge e is interpreted in the standard way, denoted e{i',j'): 

- Ri = {i',j’) e Ri 

— {Ri V R 2 ) = {i' ,3) G Riij U 

— (i?i A R 2 ) = G -^1.3 I? 

- (-Ai) = (z',/)e(G, x;7,)\i?i,^. 

An edge between i,j is satisfied by a configuration i',j' if e{i',j') is true (this 
is generalised to n-ary relations). A spatial constraint graph G = {N,E) is 
consistent or satisfiable if there exists a configuration s of N that satisfies all 
edges in E, denoted G(s); this is referred to as the consistency task. Graph G' is 
a consequence of, or implied by, G if every spatial configuration that satisfies G 
also satisfies G'. This is the sufficiency task (or entailment) that we commonly 
apply to constructive proofs, where the task is to prove that objects and relations 
in G are sufficient for ensuring that particular properties hold in G'. 

Given graph G, two key questions are (1) how to give meaning, or interpret, 
the spatial relations in G, and (2) how to efficiently determine consistency and 
produce instantiations of G. That is, we need to adopt a method for spatial 
reasoning. 

® That is, constructive methods may fail in building a consistent solution, and iterative 
root finding methods may fail to converge. 
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(a) initial configuration (b) rotation (c) translation 


(c) translation 


Fig. 4. Topological relations between four spheres maintained after various affine trans¬ 
formations. 

3.1 An Example of the Basic Concept 

A key insight is that spatial configurations form equivalence classes over qualita¬ 
tive relationships based on certain affine transformations. For example, consider 
the spatial task of determining whether five same-sized spheres can be mutu¬ 
ally touching. Suppose we are given a specific numerically defined configuration 
of four mutually touching spheres as illustrated in Figure 4(a), and we prove 
that it is impossible to add an additional mutually touching sphere to this con¬ 
figuration. That is, let si,...,S 4 be unit spheres (radius 1), centred on points 



respectively. 


According to Equation 1, Si,..., S 4 are mutually touching. We prove that a fifth 
same-sized, mutually touching sphere cannot be added to this configuration by 
determining that the corresponding system of polynomial constraints is unsatisfi- 
able (the system consists of four constraints with three free variables Xs ^, , Zs ^, 
by reapplying Equation 1 between 55 and each other sphere, e.g. si touches S 5 


is + y% + z\ 


.2 

'S5 


4). 


Now consider that we apply an affine transformation to the original configuration 
such as rotation, translation, scaling, or reflection, as illustrated in Figures 4(b) 
and 4(c). After having applied the transformation, it is still impossible to add 
a fifth mutually touching sphere, because the relevant qualitative (topological) 
relations are preserved under these transformations. Thus, when we proved that 
it was impossible to add a fifth same-sized mutually touching sphere to the 
original given configuration, in fact we proved it for a class of configurations, 
specifically, the class of configurations that can be reached by applying an affine 
transformation to the original configuration. Now, when determining consistency 
of graphs of qualitative spatial relations, we are not given any specific spatial 
configurations to work with (i.e. complete absence of numerical information), 
and instead need to prove consistency over all possible configurations. 

The key is that, each time we ground and constrain variables, we are eliminat¬ 
ing a spatial symmetry from our partially defined configuration. If we maintain 
knowledge about symmetries that certain object types have (e.g. spheres have 
complete rotational symmetry) then we can judiciously “trade” symmetries for 
unbound variables in our polynomial encoding at a purely symbolic level. Im¬ 
portantly, rather than having to compute symmetries or undertake any complex 
symmetry detection procedure, we are instead building knowledge about space 
and spatial properties of objects into the spatial solver at a declarative level. 
Thus, we are able to efficiently reason over an infinite set of possible config¬ 
urations by incrementally pruning spatial symmetries based on commonsense 
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knowledge about space, and this pruning is exploited by eliminating and con¬ 
straining variables in the underlying polynomial encoding. 


3.2 Theoretical Foundations for Symmetries 

Due to the parameterisation of objects, spatial configurations are embedded in 
n-dimensional Euclidean space K" (1 < < 3) with a fixed origin point. Let 

V,W be Euclidean spaces in M”, each with an origin. Given vectors x,y and 
constant k, a linear transformation / is a mapping V ^ W such that 

/(x + y) = f(x) + f{y) (additive) 
f(kx) = kf{x) (homogeneous) 

An affine transformation /' is a linear transformation composed with a transla¬ 
tion. It is convenient to represent a linear transformation on vector x as a left 
multiplication of a d x d real matrix Q, and translation as an addition of vector t, 
f'{x) = Qx + t. We denote a transformation T applied to a spatial configuration 
of objects s as Ts. 

We distinguish particular classes of transformations with respect to the quali¬ 
tative spatial relationships that are preserved, for example, in the following 
matrices represent rotation by 9, uniform scaling by fc > 0, and horizontal re¬ 
flection, respectively: 

/ cos{6) —sin{9)\ f k 0\ /—10\ 
sinlo) cos(9) J ’ yO k J ’ y 0 1 J ' 

Given transformation T we annotate it with its type c G C, e.g. C = {translate, 
rotate, scale, reflect} as T°. Each spatial relation R belongs to a class of relations 
in Rel, such as topology, mereology, eoincidence, relative orientation, distanee. 
Let Sym be a function Sym : Rel —>■ 2*^ that represents the classes of trans¬ 
formations that preserve a given class of spatial relations. The Sym function is 
our mechanism for building knowledge about spatial symmetries into the spatial 
reasoning system. Let Relc be the set of classes of the spatial relations that are 
used in the spatial constraint graph G, and let Syrng = H^EgReiG 

The following formal Gondition on Sym^ states that transformations (applied to 
the embedding space) define equivalence classes of configurations with respect 
to the consistency of spatial constraint graphs. When satisfied, this condition 
provides a theoretically sound foundation for symmetry pruning. 

Condition 1 Given spatial constraint graph G, configuration s, and affine trans¬ 
formation T“ with c G SymQ then G{s) is true if and only if G{T'^s). 

3.3 Polynomial Encodings for Spatial Relations 

In this section we define a range of spatial domains and spatial relations with 
the corresponding polynomial encodings. While our method is applicable to a 
wide range of 2D and 3D spatial objects and qualitative relations, for brevity 
and clarity we primarily focus on a 2D spatial domain. Our method is readily 
applicable to other 2D and 3D spatial domains and qualitative relations, for 
example, as defined in [4, 9, 10, 28, 29, 33, 34]. 


10 


Carl Schultz and Mehul Bhatt 


Helation 

Polynomial Encoding 

left of 

(point p, segment Sab) 

{Xb - Xa.){Vp - Vcl) > (Vb - ya){Xp - Xp) 

collinear 

(point p, segment Sab) 

(xb - Xa)(yp - ya) = (vb - Va)(Xp - Xp) 

right or collinear 
(point p, segment Sab) 

(Xb - Xa){yp - yp) < (yb - yp){xp - Xp) 

parallel 

(segments Sab^Scd) 

(yb - yp)(xd - Xp) = (yd - yp)(xb - Xp) 

coincident 

(point p, segment Sab) 

collinear(p, Spb) A Xp G lXp,Xb\ /\yp & [yp, yb\ 

coincident 
(point p, circle c) 

(xp - Xp)^ + (yp - yp)^ = Tp 

inside 

(point p, rectangle a) 

(0 < (p — Pla) ' < 'i^a)A 

(0 < (p — Pla) ■ '^a < ^a) 

intersects 

(point p, rectangle a) 

(0 < (p — Pla) ' '^a < ■hb'a)A 
(0 < (p — Pla) ■ '^'a ^ 

boundary 

(point p, rectangle a) 

intersects(p, a) A ^ inside(p, a) 

outside 

(point p, rectangle a) 

—1 intersects(p, a) 

concentric 
(rectangles a, b) 

hiPl’c, - Plo) +Plp = iiP3b - Pit) + Pit 

part of 

(rectangles o, b) 

/\i—i 4 intersects(pi^ , 6 ) 

proper part 
(rectangles a, b) 

^ equals(a, fe) A part of(a,fe) 

boundary part of 
(rectangles a, b) 

Ai=i ...4 boundary(pi^, b) 

discrete from 
(rectangles a, b) 

Vi=i ...4 ( right_or_collinear(a, 

V right or collinear(t>, (p^^ , p())) 

partially overlaps 
(rectangles o, b) 

3pi G ( inside(p.j, a) A inside(pi ,b))A 

3pj G ( inside(pj, a) A outside(pj , &)) A 

3pk G ( outside(pfc, a) A inside(pfc, b)) 


Table 1. Polynomial encodings of qualitative spatial relations. 


— a point is a pair of reals x, y 

— a line segment is a pair of end points pi,P 2 {pi 7 ^ P 2 ) 

— a rectangle is a point p representing the bottom left corner, a unit direction 
vector V defining the orientation of the base of the rectangle, and a real width 
and height w,h (0 < ui, 0 < h); we can refer to the vertices of the rectangle: 
let v' = {—yv, Xv) be v rotated 90° counter-clockwise, then pi = p = P5,P2 = 
wv +pi,p 3 = wv + hv' + Pi, P4 = hv' + Pi 

— a circle is a centre point p and a real radius r (0 < r) 


We consider the following spatial relations: 

Relative Orientation. Left, right, collinear orientation relations between points 
and segments, and parallel, perpendicular relations between segments [23]. 

Coincidence. Intersection between a point and a line, and a point and the bound¬ 
ary of a circle. Also whether the point is in the interior, outside or on the bound¬ 
ary of a region. 

Mereology. Part-whole relations between regions [37]. 


Table 1 presents the corresponding polynomial encodings. Given three real vari¬ 
ables v,i,j, let: 

V G [i,j] = i<v<j\/j<v<i. 
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h 

h. 


'h 

(d) rotation 


(a) configuration of (b) vertical reflection 
points and lines 


(c) x-scallng 


Fig. 5 . Affine transformations preserve point coincidence, parallelism, and ratios of 
distances along parallel lines. 

Determining whether a point is inside a rectangle is based on vector projection. 
Point p is projected onto vector v by taking the dot product. 


i^^p^yp) * i^^v^Vv) — 3^p^v “t VpVv 


Given point a and rectangle b, we translate the point such that the bottom left 
corner of b is at the origin, project Pa on the base and side vectors of b, and 
check whether the projection lies within the width and height of the rectangle, 

0 < {pa -pij -Vb < Wb 

0 < (Pa -Pu) -Vb < hb 

Convex regions a, b are disconnected iff there is a hyperplane of separation, i.e. 
there exists a line I such that a and b lie on different sides of 1. This is the basis 
for determining the discrete from relation between rectangles. 

3.4 Formalising Knowledge about Symmetries 

In this section we formally determine the qualitative spatial relations that are 
preserved by various affine transformations. A fundamental property of affine 
transformations is that they preserve (a) point coincidence (e.g. line intersec¬ 
tions), (b) parallelism between straight lines, and (c) proportions of distances 
between points on parallel lines [26]. For example, consider the configuration of 
points Pa,Pb,Pc,Pd and lines li,l 2 ,h in Figure 5: (a) we cannot introduce new 
points of coincidence between lines by applying transformations such as transla¬ 
tion, scaling, reflection, and rotation. Conversely, if two lines intersect, then they 
will still intersect after these transformations; (b) lines li,l 2 are parallel before 
and after the transformations; lines li,l 3 are always non-parallel; (c) the ratio of 
distances between collinear points Pa,Pb,Pc is maintained; formally, let Sy be the 
segment between points Pi,Pj and let |sijj be the length of the segment. Then 
the ratio in Figure 5 is the same before and after the transformations. 

Based on these properties we can determine the transformations that preserve 
various qualitative spatial relations. 

The properties of affine transformations and the geometric objects that they preserve 
are well understood; further information is readily available in introductory texts 
such as [26]. Our key contribution is formalising and exploiting this spatial knowledge 
as modular and extensible common-sense rules in intelligent knowledge-based spatial 
assistance systems. 
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Theorem 1 The following qualitative spatial relations are preserved under trans¬ 
lation, scale, rotation, and reflection (applied to the embedding space): topology, 
mereology, coincidence, collinearity, line parallelism. 

Proof. By definition, affine transformations preserve parallelism with respect to 
qualitative line orientation, and point coincidence. Due to preservation of point 
coincidence and proportions of collinear distances by affine transformations, it 
follows that mereological part of and topological contact relations between re¬ 
gions are preserved, i.e. if a mereological or topological relation changes between 
regions a, b, then by definition there exists a point p coincident with a such 
that the coincidence relation between p and b has changed; but this cannot oc¬ 
cur as point coincidence is maintained with affine transformations by definition, 
therefore mereological and topological relations are also maintained. 

The interaction between spatial relations and transformations is richer than we 
have space to elaborate on here, i.e. not all qualitative spatial relations are 
preserved under all affine transformations; orientation is not preserved under 
reflection (e.g. Figure 5(b) gives a counter example), distance is not preserved 
under non-uniform scaling. To summarise, we formalise the following knowledge 
as modular commonsense rules in CLP(QS): point-coincidence, line parallelism, 
topological and mereological relations are preserved with all affine transforma¬ 
tions. Relative orientation changes with reflection, and qualitative distances and 
perpendicularity change with non-uniform scaling. Spheres, circles, and rect¬ 
angles are not preserved with non-uniform scaling, with the exception of axis- 
aligned bounding boxes. 


“Trading” Transformations Symmetries are used to eliminate object vari¬ 
ables. As a metaphor, unbound variables are replaced by constant values in 
“exchange” for transformations. We start with a set of transformations that can 
be applied to a configuration: translation, scaling, arbitrary rotation, and hor¬ 
izontal and vertical reflection. We can then “trade” each transformation for an 
elimination of variables. Each transformation can only be “spent” once. Theo¬ 
rem 2 presents an instance of such a pruning case. 

Table 2 presents a variety of different pruning cases for position variables and 
the associated combination of transformations, as illustrated in Figure 6 .^^ Some 
cases require more than one distinct set of parameter restrictions to cover the 
set of all position variables due to point coincidence being preserved by affine 
transformations. For example, consider case (f): all pairs of points pi,p 2 can 
be transformed into any other pair of points Pi,Pj by translation, rotation, and 
scaling, iff pi = P 2 Pi = Pj- Thus, to cover all pairs of possible points, we need 
to consider two distinct parameter restrictions: pi = pj and pi Pj', we refer to 
these as subcases. 

Many further pruning cases are identifiable. For example, a version of case (i) 
can be defined without reflection by requiring more sub-cases where C 4 > C 2 and 
C 4 < C 2 . Case (i) can be extended so that all six coordinates of three points 
are grounded if we also “exchange” the skew transformation (e.g. applicable to 
object domains like triangles or points). 
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All cases have been verified using Reduce as presented in Theorem 2. 
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Theorem 2 Any pair of object position variables (a;i,t/i), (a;2, j/2) can be trans¬ 
formed into any given position constants (ci,C2), (03,C2) such that (ci =030 
(xi = X 2 /\ yl = y2)) by applying: an xy-translation, a rotation about the origin 
in the range (0, 27r), and an x-scale. 

Proof. The corresponding expression has been verified using the Reduce system 
(Redlog quantifier elimination) [15]; all variables are quantified over reals. 


'iciic2ic3ixiiy\ix2iy2 

(ci = C 3 -H- (*1 = *2 A t/1 = 1/2)) -f-/- 

(0 < Sx) A (d^ + dy = 1) A 



We can use this pruning case on any spatial constraint graph G where the graph’s 
spatial relations are preserved by translation, rotation, and scaling. Given graph 
G = {N,E), the following Algorithm applies the pruning case, with selected 
constants ci = 0 , C 2 = 0 , and C 3 = 1 or C 3 = 0 : 

1. select object position variables Pi,P 2 from nodes in N 

2. copy G to create Gi, G 2 

3. in Gi set pi = (0, 0 ),p 2 = (1, 0) (case ci ^^ 03 ) 

4. in G 2 set pi = (0,0),P 2 = ( 0 , 0 ) (case Ci = 03 ) 

5. if the task is: 

(a) consistency of G then solve V^=i Gi(s) 

(b) sufficiency, G —)■ G', then solve 
ALi-3s(G.(s)A-G'(s)) 

In Step 1 any pair of objects can be selected for which their position variables 
will be grounded; we also employ policies that target computationally costly 
subgraphs (for example, pairs of non-equal circles that share a boundary point 
are often good candidates for this pruning case). Having eliminated free variables 
from the system of polynomial constraints, the constraints are significantly more 
simple to solve. Due to the double exponential complexity 0{ci ^) reducing n 
has a significant impact on performance; the system may even collapse from 
nonlinear constraints to linear (solvable in 0 (c”)) or constants. 


3.5 Combining Symmetry Pruning with Graph Decomposition 

In certain cases, spatial constraint graphs can be decomposed into subgraphs 
that can be solved independently. For example, subgraphs Gi,G 2 can be inde¬ 
pendently solved if all objects in subgraph Gi are either: 

— disconnected from all objects in subgraph G 2 ; 

— a proper part of some object in G 2 ; 

— left of some segment in G 2 ; 

— only related by relative size to some object in G 2 , and so on. 
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X}=C} 


(a) 


xi=ci \ i 

yi=c2 \ i 

o 

•Y i 

Xl=Cl X2=C2 

(b) (c) 




II 

X2 

=C3 




Xi= 

Cl 

' X1=C1 

yi=C 2 


X 

Xl=Cl 

yi=C2 


(d) (e) (f) 





y3=C3 


y3=C4 

y3=C4 








X} = 

Cl 

X2 

X 

X1=C1 

=C2 yi=C2 

X2 

^ Xl=ci 

=C3 yi=C2 

X2=C3 ■ 

y2=C2 


(g) (h) (i) 


Fig. 6. Cases for pruning position parameters. 

In such cases we can reapply spatial symmetry pruning in each independent 
sub-graph; this commonsense spatial knowledge is modularly formalised within 
CLP(QS). For example, consider Proposition 22 of Book I of Euclid’s Elements 
(Figure 7): 


Constructing a triangle from three segments. Given three line segments lat, led, 
/e/, draw a line through four collinear points pi,... ,p 4 such that |(pi,p 2 )| = |Gh|j 
l(P 2 ,P 3 )| = \lcd\, \iP 3 ,P 4 :)\ = \lef\- Draw circle Ca centred on p 2 , coincident with 
Pi. Draw circle Cb centred on pa coincident with p 4 . Draw ps coincident with 
Co and Cb- The triangle P 2 ,P 3 ,P 5 has side lengths such that |(p 2 ,P 3 )| = |^cd|j 

l(P3,P5)| = \lef\, \{P5,P2)\ = \lab\- 


In this example, the three segments lab, led, lef and the remaining objects are 
only related by the distances between their end points. That is, the relative 
position and orientation of lab, led, lef is not relevant to the consistency of the 
spatial graph; we only need to explore all combinations of segment lengths. Thus 
the solver decomposes the graph into four sub-graphs: (1) lab (2) led (3) lef, and 
(4) pi,... ,p 5 , Co, Cfe. In subgraphs (I),(2),(3) it “trades” translation and rotation 
to ground Pa = Pe = Pe = (0,0), and Pb = Vd = Vf = 0 and keeps a;-scale to cover 
all possible combinations of segment lengths, i.e. Xb,Xd,Xf are free variables. In 
subgraph (4) CLP(QS) applies the pruning case of Theorem 2 by grounding 
Pi, Pi- 
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Case 

Parameter restrictions 

'iVaded transformations 

a 

Xi — Cl 

fc-translate 

b 

— Ci.y-i — C2 

fcy-translate 

c 

X\ — Ci,X2 — C2, (i) Cl 7 ^ C2 (ii) Ci = C2 

fc-translate, rotate tt, T-scale 

d 

— Cl, y2 — C2 

fcy-translate 

e 

— Ci,yi — C2,X2 — C3, 

(i) Cl ^ C3 (ii) Cl = C3 

aiy-translate, rotate tt, rc-scale 

f 

— Cl, yi — C2,X2 — C3,y2 — C 2 , 

(i) Cl 7 ^ C3 (ii) Cl = C3 

fcy-translate, rotate ( 0 , 27 r), x-scale 

g 

Xi = Ci,X2 — C2,y3 — C 3 
(i) Cl ^ C 2 (ii) Cl = C 2 

xy-translate, rotate tt, x-scale 

h * 

xi — Cl, yi — C2,X2 — cz,y3 — C4 

(i) Cl 7 ^ C3 A C2 7 ^ C4 (ii) ci — /\ c2 ^ C4 

(iii) Cl 7^ C3 A C2 = C4 (iv) ci = C3 A c 2 = C4 

xy-translate, rotate tt, xy-scale, 
2/-reflect 

i * 

Xl = Ci,yi — C2,X2 = C3,y2 — C2,y3 — C4 
(i) Cl 7 ^ C3 A C2 7 ^ C4 (ii) ci / C3 A C2 = C4 
(iii) Cl = C3 A C2 = C4 

xy-translate, rotate ( 0 , 27 r), xy-scale, 
y-reflect 


Table 2. Cases for pruning parameters for one position point (a,b), two position points 
(c-f), three position points (g-i)- Cases marked with * require arbitrary scaling (i.e. both 
uniform and non-uniform). 


?- Sab = segment (Pa,Pb) , Scd = segment (Pc , Pd) , Sef = segment (Pe,Pf ), 
I S12 = segment(PI,P2), S23 = segment(P2,P3), S34 = segment(P3,P4), 
I S14 = segment(PI,P4), 

I coincident(P2,S14), coincident(P3,S14), 

I 

I length_dim(equal, S12, Sab), 

I length_dim(equal, S23, Scd), 

I length_dim(equal, S34, Sef), 

I 

I Ca=circle(P2,_), coincident(PI,Ca), 

I Cb=circle(P3,_), coincident(P4,Cb), 

I coincident(P5,Ca), coincident(P5,Cb), 

I S35 = segment(P3,P5), S52 = segment(P5,P2), 

I 

I y,X triangle 

I length_dim(equal, S23, Scd), 

I length_dim(equal, S35, Sef), 

I length_dim(equal, S52, Sab), 
true . 

I (length_dim(not_equal, S23, Scd); 

I length_dim(not_equal, S35, Sef); 

I length_dim(not_equal, S52, Sab)), 
false . 



Fig. 7. Constructing a triangle by decomposing the spatial constraint graph. 


4 Application-Driven Use Cases 


We present problem instances in the classes of mereology, ruler and compass, 
and contact. Table 3 presents the experiment time results of CLP(QS) using 
symmetry pruning compared with existing systems: z3 SMT solver, Redlog real 
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quantifier elimination (in the Reduce computer algebra system) [15], and the re¬ 
lation algebraic qualitative spatial reasoners GQR [16] and SparQ [39] . CLP(QS) 
uses z3 to solve polynomial constraints (after our pruning), thus z3 is the most 
direct comparison. Experiments were run on a MacBookPro, OS X 10.8.5, 2.6 
GHz Intel Gore i7. The empirical results show that no other spatial reason¬ 
ing system exists (to the best of our knowledge) that can solve the range of 
problems presented in this section, and in cases where solvers are applicable, 
GLP(QS) with spatial pruning solves those problems significantly faster than 
other systems. 


Problem 

CLP(QS) 

z3 

Redlog 

GQR 

SparQ 

Aligned Concentric 

6.831 

47.651 

time out 

n/a 

n/a 

Boundary Concentric 

2.036 

time out 

time out 

n/a 

n/a 

Mereologically Concentric 

0.105 

0.373 

time out 

n/a 

n/a 

Angle Bisector 

0.931 

time out 

time out 

n/a 

n/a 

Sphere Contact 

0.004 

time out 

time out 

fail 

fail 


Table 3. Time (in seconds) to solve benchmark problems using CLP(QS) with prun¬ 
ing compared to z3 SMT solver, Redlog (Reduce) guantifier elimination, and GQR and 
SparQ relation algebraic solvers. Time out was issued after a running time of 10 min¬ 
utes. Failure (i&\\) indicates that the incorrect result was given. Not applicable (nfa.) 
indicates that the problem could not be expressed using the given system. 


4.1 Spatial Theorem Proving: Geometry of Solids 

Tarski [36] shows that a geometric point can be defined by a language of mereo- 
logical relations over spheres. The idea is to distinguish when spheres are concen¬ 
tric, and to define a geometric point as the point of convergence. Borgo [8] shows 
that this can be accomplished with a language of mereology over hypercubes. We 
will use GLP(QS) to prove that the definitions are sound for rotatable squares. 

As a preliminary we need to determine whether the intersection of two squares 
is non-square (Figure 8) [34] . Given two squares a,b, the intersection is non¬ 
square if a partially overlaps b (Table 1) and either (a) a and b are not aligned, 
^ or (b) the width and height of the intersection are not equal, wj ^ hi, 
such that 

wi = min(w • p 2 „,• P 2 J - max(u • u • pij 

hi = min(u' • p4„, u' • p 4 J - max(u' • , u' • J 


Aligned Concentric. Two squares A, B are aligned and concentric if: A is part of 
B and there does not exist a square P such that (a) P is covertex with B, and 
(b) the intersection of P and A is not a square (Figure 9). 

I aligned_concentric(A,B) 

I mereology (part_of , A,B), 

I not(( 

I P=square , 

I covertex(P,B) , 

I intersection(non_square,A,P) . 

I )). 
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B 



I 



A 





. 


Wi 

(b) hi 7 ^ wi 


Fig. 8. Intersection I of squares A, B is non-square. 



(a) A and B are 
non-concentric 



(b) A and B are 
concentric 



(c) A is bound¬ 
ary concentric 
with B. 



(d) A and B are 
rotated concen¬ 
tric. 


Fig. 9. Characterising aligned (a),(b) and rotated (c),(d) concentric squares using 
mereology (reproduced from [Borgo, 2013]) 

We use CLP(QS) to prove that the definition is sufficient, by contradiction; two 
squares are covertex if they are aligned and share a vertex, and the relation 
concentric is the geometric definition of concentricity in Table 1 that is used to 
evaluate the mereological definition of concentricity: 

?- A=squ8ire , 

I B=squ£a:e , 

I allgrLed(A,B) , 

I not_concentric(A,B) , 

I aligned. concentricCA.B) . 
false . 

Boundary Concentric. Square A is boundary concentric with square B if: A is 
proper part of B and there does not exist a square Z such that (a) Z is proper 
part of B (b) A is part of Z, and (c) Z is not part of A (Figure 9). 


I boundary_concentric(A,B) 

I mereology (proper.part , A , B) , 

I not(( 

I Z = square 

I mereology (proper.part , Z,B), 

I mereology (part.of , A,Z), 

I mereology (not_part_of , Z,A), 

I )). 


?- A=square , 

I B=square , 

I not_aligned(A,B) , 

I not_concentric(A,B) , 

I boundary_concentric(A,B) . 
false. 
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Mereologically Concentric. Squares A, B are mereologically concentric if: A, B 
are aligned concentric or there exists Q such that (a) Q is boundary concentric 
with B and Q is aligned concentric with A or (b) Q is boundary concentric with 
A and Q is aligned concentric with B. 

Having proved the mereological definitions of aligned and boundary concentricity, 
we can replace these with more efficient geometric definitions from Table 1 when 
proving mereological concentricity. 

I g_aligned_concentric(A,B) 

I aligned (A,B) , concentric(A,B) . 

I 

I g_boundary_concentric(A,B) 

I mereology (boundary_part , A,B). 

I m_concentric(A,B) 

I g_aligned_concentric(A,B) 

I ; 

I Q = square(_,_,_) , 

I g_boundary_concentric(Q,B) , 

I g_aligned_concGntric(Q,A) 

I ; 

I Q = squarG(_,_,_) , 

I g_boundary_concentric(Q,A) , 

I g_aligned_concGntric(Q,B) . 

?- A=squarG , 

I B=squarG , 

I not_concentric(A,B) , 

I m_concentric(A,B) . 
false. 


4.2 Didactics: Ruler—Compass and Contact Problems 


Angle Bisector. Let la,lb be line segments that share an endpoint at p. Draw 
circle c at p. Circle c intersects la at Pa and 4 at pb. Draw circles Ca at pa and 
circle Cb at pb such that p is coincident with both Ca,Cb. Circles Ca and Cb intersect 
at p and p^. The line segment from p to p^ bisects the angle between la and 4 
(Figure 10). 

We use CLP(QS) to prove that the definition is sufficient. The relation bisects is 
used for evaluation by checking if the midpoint of {pa,Pb) is collinear with 4 (i-e. 
idealised rulers cannot directly measure the midpoint of a line). An interactive 
diagram is then automatically generated that encodes the specified program 
(using the FreeCAD system); see Figure 11. 

?- La = segment (P, _), Lb = segment (P , 

I orientation(not_parallel, La, Lb), 

I C = circle(P,_), 

I coincident(Pa,La),coincident(Pa,C), 

I coincident(Pb,Lb),coincident(Pb,C), 

I Ca=circle(Pa,_), Cb=circle(Pb,_), 

I coincident(Pa,Ca),coincident(Pa,Cb), 

I coincident(Pc,Ca),coincident(Pc,Cb), 

I not_equal(P,Pc), 

I Lc=segment(P,Pc), 

I bisects(segment(Pa,Pb),Lc). 
true . 

I not_bisects(segment(Pa,Pb),Lc). 
false . 
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Fig. 10. Ruler and compass method for angle bisection. Line Lc bisects the angle 
between lines La, Lb- 



(a) (b) (c) 


Fig. 11. Interactive diagram encoding the student’s constructive proof of Euclid’s an¬ 
gle bisector theorem; as the student manipulates figures in the diagram, the other 
geometries are automatically updated to maintain the specified qualitative constraints. 

Sphere contact. Determine the maximum number of same-sized mutually touch¬ 
ing spheres (Figure 4 - note that no numeric information about the spheres is 
given in this benchmark problem). 

?- spJiereJJs^ (5, Spheres), 

I size (equal, ^roup (Spheres), 

I ropcao^y(rcc8(ec), ^roup (Spheres) ). 
raise. 


5 Conclusions 


Afhne transformations provide an effective and interesting class of symmetries 
that can be used for pruning across a range of qualitative spatial relations. 
To summarise, we formalise the following knowledge as modular commonsense 
rules in CLP(QS): point-coincidence, line parallelism, topological and mereolog- 
ical relations are preserved with all affine transformations. Relative orientation 
changes with reflection, and qualitative distances and perpendicularity change 
with non-uniform scaling. Spheres, circles, and rectangles are not preserved with 
non-uniform scaling, with the exception of axis-aligned bounding boxes. Our al¬ 
gorithm is simple to implement, and is easily extended to handle more pruning 
cases. 

Theoretical and empirical results show that our method of pruning yields an 
improvement in performance by orders of magnitude over standard polynomial 
encodings without loss of soundness, thus increasing the horizon of spatial prob¬ 
lems solvable with any polynomial constraint solver. Furthermore, the declara- 
tively formalised knowledge about pruning strategies is available to be utilised in 
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a modular manner within other knowledge representation and reasoning frame¬ 
works that rely on specialised SMT solvers etc, e.g., in the manner demonstrated 
in ASPMT(QS) [38], which is a specialised non-monotonic spatial reasoning sys¬ 
tem built on top of answer set programming modulo theories. 
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